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Abstract 

The effect of the inclusion of ground state correlations into the QRPA equation of 
motion for the two-neutrino double beta (/3/3 2v ) decay is carefully analyzed. The resulting 
model, called renormalized QRPA (RQRPA), does not collapse near the physical value of 
the nuclear force strength in the particle-particle channel, as happens with the ordinary 
QRPA. Still, the Pfav transition amplitude is only slightly less sensitive on this parameter 
in the RQRPA than that in the plain QRPA. It is argued that this fact reveals once more 
that the characteristic behaviour of the Pfav transition amplitude within the QRPA is not 
an artifact of the model, but a consequence of the partial restoration of the spin-isospin 
SU(4) symmetry. It is shown that the price paid for bypassing the collapse in the RQRPA 
is the violation of the Ikeda sum rule. 
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1 Introduction 



The quasiparticle random phase approximation (QRPA) is the most frequently used nu- 
clear structure method for evaluating (3/3 rates both for the two-neutrino decay mode 
{PPiv) and for the neutrinoless mode (J3(3o v )- The general feature of this method is that 
the resulting nuclear matrix elements M.2u and M.q v turn out to be highly sensitive to 
the particle-particle force in the S — 1, T — channel [|l], |, || ||, ||, ||]. Furthermore, 
the QRPA collapses very close to the physical value for this force. One thus may suspect 
that this method yields relatively small values of M.i v simply because the approximation 
breaks up. In other words, the smallness of M.% v in the QRPA could be just an artifact 
of the model. 

Several modifications of the QRPA have been proposed in order to amend the above 
behavior in a qualitative way, including higher order RPA corrections J7|, nuclear de- 
formation 0, single-particle self-energy BCS terms ||, particle number projection [|IIJ] 



and the proton-neutron pairing (]]J. Yet, none of these inhibits the collapse, which is 
just the famous "phase transition" where the RPA develops zero or imaginary frequency 
solutions. One should remember in this connection that in the derivation of the ordinary 
RPA it is assumed that we can replace, when evaluating the equations of motion, the 
RPA (correlated) ground state by the Hartree-Fock ground state. Discrepancies due to 
this replacement will obviously get more serious the more significant are the ground state 
correlations. Attempts to correct this have been made by embodying the effect of GSC 
in the RPA equations of motion [|l^, [13|, [L4|, [15|]. The corresponding formalism, named 
self-consistent or renormalized RPA (RRPA), generally yields better results in the sense 
that the instability develops at a larger interaction strength or is fully taken away. 

Recently, the renormalized QRPA (RQRPA) has been applied to the fifiiv decays by 
Toivanen and Suhonen [|16| (in 100 Mo) and by Schwieger, Simkovic and Faessler |17j (in 
76 Ge, 82 Se, 128 Te and 130 Te nuclei). We have also discussed the 100 Mo 2z/-decay but only 
in the framework of a schematic model ||18|| . In this work we present a detailed study of 
several f3f3 decaying nuclei within the RQRPA. Particular attention is given to the non- 
conservation of the non energy weighted sum rule (called Ikeda sum rule for J n = 1 + ) 
within this approximation. The self-consistency between the residual interaction and the 
mean field || is also addressed. 
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2 Formalism 



The RQRPA formalism has been presented in several articles ITR jTR IT7L ITS] and we will 



therefore only sketch the main formulae needed in the discussion. 

2.1 RQRPA Equation and the Transition Moments 

The excited states \XJM) are built by the action of the charge-exchange operators 

Q\\JM) = £ [x pn (XJ)Al n (JM) - Y pn (\J)A pn (JM)] , (1) 

pn 

on the QRPA correlated ground state |0). Here A\ n (JM) = [aJc^] JM , A pn (7M) = 
(— 1) J+M A pn { J — M), and a\ are quasiparticle creation operators for protons (t = p) and 
neutrons (t = n). The amplitudes X and Y and the eigenvalues u)\ satisfy the RQRPA 
equations 



A(J) B(J) \ / X(AJ) \_ ( X(AJ) \ m 
B*(J) A*(J) i I Y(AJ) J AJ I -Y(AJ) i ' 1 ' 



where 



X pn (XJ) = X pn (XJ)DH 2 and Y pn (AJ) = F pn (AJ)^/ 2 (3) 
are the renormalized amplitudes, 

-Dpn — l—Mp- H n (4) 

and 



M = j^(Q\[a\a,m (5) 

are the quasiparticle occupations (J = V2J + 1). Here we are assuming that the single 
quasi-particle density is diagonal, (0| [ajap] ^) = 5«'(0|[a|o!f] |0). The submatrices A(J) 
and B(J) are 
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hpn, P >n'{J) = (e p + e n )<W<W + D^ n 2 [F(pn, p'n', J) (u p v n u p >v n , + v p u n v pl u n ) 

+ G(pn,p'n\ J)(u p u n u p iu n i + v p v n v p >v n >)\ D^,, 
Bpn, P 'n'{J) = D p n [F(pn,p'n, J)(v p u n u p ,v n > + u p v n v p >u n >) 

1 /2 

- G(pn,pn, J)(u p u n v p >v n > + v p v n u p >u n >)} DJ n ,, 



(6) 



where F and G are the usual particle-hole (ph) and particle-particle (pp) coupled two-body 
matrix elements. 

The QRPA ground state is determined by the condition 



tt{\JM)\0) = 0. 
Writing the ground state in the form 

\0) = N e s \BCS), 
using equation fl7|) and making the quasi-boson approximation, we find 



S — - ^2 J 1 D p n^ 2 C pnp ' n i(J)D p , n , A pn (J)A p , n ,(J) 

pnp'n'J 

where the matrix C is the solution of 



(7) 



(9) 



^X; n (AJ)C^ n ,(J) = Y;, n ,(AJ), for all XJ. 

pn 

From these equations one gets that 

K = Y.J 2 h 2 Dpn'\ypn'{\J)\\ 

\Jn> 

K = Y,j 2 3n 2 Dpi n \YM\J)\ 2 . 

XJp' 

To evaluate the transition matrix elements 



(aj||o ± (j)||o) = (o| n(\j),o ± (j) |o), 



o ,~, 



the P ± decay operators 



(10) 



:n) 



(12) 
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± (JM) = Y J 0{JM;i)t ± (i); (with t-\n) = \p)) 

i 

are expressed in the form 



where 



with 



Of(JM) 
Ol\JM) 



± {JM) = Of(JM) + Of{JM), 



pn 

£ [^ pn {J)Bl n {JM) + (-1)^(7)^(7^) 



(13) 



(14) 



(15) 



Kn(J) = -(-V J j~\u n (p\\0(J)\\ny, A pn (J) = -J- 1 u p v n (p\\0(J)\\n), 
A+(J) = -(-l) J J^Vn(p||0(J)||n)*, k-(J) = J- l u p u n l$\\0{J)\\n), (16) 



Bl n (JM) = [ala n ] JM and B pn (JM) = (-l) J+M B pn (J - M). Notice that, whereas the 
operator O^(JM) creates and destroys a np quasi-particle pair, the scattering operator 
O^(JM) conserves the number of quasi-particles. 



The evaluation of (12) goes now via the following relationships: 



5pp* S nn i D pn , 



r\V\[A pn {J),Al, n ,{J)\ |0) 
J- x (0| \B pn {J), Bl, n ,{J)]° |0> = 8 pp ,8 nn ,{U n - U p ), 



J-^OI k„(J), Bl, n ,(J)]° |0> = -8 nn ,j^(0\ 



J~ L {0\ AUJ),B p , n ,{J) |0> = -VJn^Ol 



o ,~, 



t t 

(X i 



|o> = o, 
° |5>^o. 



The expectation values (0| 



o ~ 



(17) 



|0) are identically zero since |0) is a superposition of 



states with equal number of neutron and proton quasi-particles (see eqs. (ffi) and (|1)). 
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Thus, only the first term in eq. (|I4| ) contributes to the transition matrix element and one 

gets 



(aj||0 ± (j)||o) = (o\[n(\j),of(j)} |o) 

= J^[A p ± n (J)X;„(AJ) + (-l) J A p ± :(^)Y; n (AJ)] J D;/ 2 (If 



pit 

The corresponding RQRPA total strengths are 



S ± (J) = J- 2 £|(AJ||O ± (J)||0>| 2 , (19) 

A 



and from (ITBD we obtain that 



S-(J)-S + (J) = T,(\KnU)\ 2 -\Kn(J)\ 2 )D P n 
pn 

= J- 2 Y.«-vl)Dp n \(p\\0{J)\\n)\\ (20) 

pn 

For the evaluation of the fifov matrix elements, instead of the two- vacua QRPA method 



introduced in ref. |19 l, we will here simply use the expression 



but solving the gap equations for the intermediate nucleus |[20|| . We have tested numeri- 
cally that both methods yield almost identical results. 

2.2 Sum Rule 

It is well known that the total (3^ strengths S±(J) 

s ± (J) = J- 2 £IM|o ± (J)||o>| 2 (22) 

V 

can be expressed in the form 

S ± (J) = (-) J J- 1 (0\[O,(J)O ± (J)}°\0}, (23) 
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when \uJ) is the complete set of excited states that can be reached by operating with 
Q T (J) on the initial state |0). It follows at once that 



S_(J) - S + (J) 



(24) 



i.e., the difference between the exact total S_(J) and S + (J) strengths equals the expecta- 
tion value of the operator (— ) J J _1 [(9 + (J), O-(J)] in the ground state |0). In particular, 
for the observables of interest here, i.e., the Fermi (F) and Gamow- Teller (GT) operators 
t± and t±a, one gets [f21|, ^ 



S_(r = 0+ 1+) - S+(J* = + , 1+) 
Moreover, if the ground state is chosen so that 



(0\N-Z\0). 



(25) 



(0\N-Z\0) = N-Z, 



one obtains the well known result 



(26) 



s_(.r = o+, i + ) - s + {r = o+, i + ) = n- z. 



(27) 



It can now be easily shown that the relations ( p3|) and ( 24]) are not valid within the 
RQRPA, and this leads to the violation of the Ikeda sum rule. In fact, using eqs. flTT| ) 
the right hand side of eq. (EM) evaluates in the RQRPA to 



(-) J J- 1 (0|[O + (J),O_(J)]°|0) = (-) J J- 1 (0|[Of(J),O 20 (J)]°|0) 

+ (-) J J- 1 (0|[O} 1 (J),O n (J)]°|0). (28) 

Calculating the first term on the right hand side of eq. (p8|) one finds that it reproduces 
the difference between the total RQRPA strengths S-(J) and S+(J) given in eq. (|20|) , 
that is 



(-) J J- 1 (0|[Of(J),O 20 (J)]°|0) = J- 2 ^n-vl)D P n\(p\\0(mn)\ 2 = S4J)-S + (J). 

pn 

(29) 
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As for the second term on the right hand side of eq. (|28| ) one obtains 



H J J- 1 (0\[01\J),0 1 _\J)} \0) = J- 2 Y.\(v\\0{J)\\n)\\ul-vl){N n - N p ). (30) 

pn 

This quantity is different from zero because the states 0^{JM)\$) are both non-null and 
orthogonal to the RQRPA model space (spanned by the states flt (AJM)|0)). 

Using the above results we can express S^(J) — S + (J) for F and GT transitions as 



5_(j^ = o + ,i + ) - s + (r = o + ,i+) = (o\n-z\o) 

J- 2 Y. MP{r = o + , i + )HI 2 (^ - v 2 p )(K-K), (si) 



pn 



and the violation of the corresponding sum rules within the RQRPA is associated with a 
nonvanishing value of the quantity 



A5(j" = o+ i + ) = (n — z) — s_{r = o+, i + ) - s+ir = o + , r 



(32) 



As discussed in Sect. 3 below, this depends quantitatively on the way the gap equations 
are solved. In fact, the usual constraints [1IBI IT?! 



N = (BCS\N\BCS) 
Z = (BCS\Z\BCS) 



2 2 
n u ni 



E"2 2 
3 r v. 



p^pi 



(33) 



lead to results different from those of ref. ITS 



AT = <0|iV|0> = Y,jX + «-<Wn], 
n 

Z=(0\Z\0) = Y,J>l + {ul-vl)N p \. (34) 
v 

In ordinary QRPA, where |0) -> \BCS), the states 0£( JM)\BCS) are null-vectors and 
the F and GT sum rules are fulfilled, i. e., the quantity corresponding in this approximation 
to eq. ( |32"D vanishes. 
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2.3 Self-consistency between the Residual Interaction and the 
Mean Field 



It is well known that in the limit of exact isospin symmetry all the S-(J n = + ) strength 
is concentrated in the isobaric analog state (IAS), there is no (3 + strength and the fifiiv 
decay is forbidden, i.e., 

s IAS = s_(r = o + ), s + (r = o + ) = o, M 2 „{r = o+) = o. (35) 

The question of which conditions have to be fulfilled in order for these relations to hold 
in the QRPA has been discussed in ref. ||. Within the RQRPA they read 

= Dli 2 A; n (r = o+)/VN^ = vW^Dli 2 j pUp v n , 

Y p f = -Dli 2 A; n (r = + )/VN^=-VN^Dli%u n v p , (36) 



as can be checked from eqs. fll6"D and fll8|). Putting these expressions into the RQRPA 
equations (0) we get 



, _ TTIAS , V P\IAS, U n*IAS 

e p + e n -u 1AS - -U jp=jn + — A p + — A n , 

u p u n 

e P + e n + u IAS = U£° +^A- s +^Ar. (37) 



Jp—Jn „. P 

V p Ui n 



where 



K AS = -\Y.3p'h l G(pp,v'v', r = + )u p ,v p ,D p , p ,, 
p' 

K AS = -^Y,Jn'Jn 1 G(rin,n , n , ,.r = + )u n ,v n ,D n , n ,, (3* 

n' 

U L= jn = E j/j^Fipn^p'n^r = 0+)(v 2 n , -v 2 pl )D p , n , 

jp>=jn> 



In deriving the eqs. (p7|) the relation 2G(pn,p'n' , J n = + ) = G(pp,p'p' , J w = + ) = 
G(nn,n'n' , J 71 = + ) has been used. We have also assumed that J 71 = + is the only 
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important multipole affecting the solution of eq. ([11]). This leads to the identity D pn 
Dpp D nn . 

Summing and subtracting the two equations in fl3"7|) we also obtain 



A IAS AIAS 

€p _|_ €n = — E 1 Vl — ? 

/\ias A IAS 

tow = U^ n + {ep-\ p )^-{e n -X n )^. (39) 

As et = A t /2vtUt, the first relation in (|39| ) will be an identity only if At = A' t AS . This 
means that in solving the gap equations the substitution 



G(tt, ff, r = 0+) -> G{tt, ff, r = 0+) Dm (40) 

has to be done. 

Also from the last relation in (|39"D we obtain 

w /AS + A p -A n = e p -e n + C^ 1 i in . (41) 

The left hand side in this equation is just the excitation energy of the IAS relative to the 
ground state of the initial nucleus, i. e., E IAS — E^. [] Further e p — e n = A c — Uj y 2j n , where 
Ac and U^j n are, respectively, the Coulomb displacement energy and the symmetry 
energy Thus, E IAS - E t will be equal to A c only if Ug? jn = Uj^ jn . But, in the BCS 
approximation, the symmetry energy reads 



U; y p =u = E vj p - 1 F(pn,pV;0 + )(^,-^ ) 

= \ E E(-WW; JT){vl - v],), (42) 
i' P =i' n JT 

and is equal to 

1 It should be remembered that in the BCS approximation the energy difference between the ground 
states of an even-even (N,Z) nucleus and an odd-odd (N — 1, Z + 1) nucleus is E N _ liZ+1 — E N Z — 
Ap + A n + Ap — A n ; also E N _ 2 ,z+2 — E NZ — 2(Ap A n ). 
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i' j, - /'„• (43) 

where 

Vjn = ^jn l jn>v 2 n ,F{nn, n'n'; + ) + ^ j" 1 j p >v 2 p ,F(nn } p'p' '; 0+) 

n' p' 

= JEE JT 2 [G(jfjf; JT = l)(2^ + v\) + G(jj'jf; JT = 0))v}], 
fo P = Yl 3p l 3p'V 2 p >F(pp, pp'; 0+) + 3p 1 3n'V 2 n ,F(pp, n'n'] + ) 

p 1 n' 

= ^EE J 2 J- 2 [G(jfjf; JT = l)(2v 2 , + v 2 ,) + G(jfjf; JT = 0))v 2 ,]. (44) 

are, respectively, the neutron and proton self-energies. This suggests that, when the 
self-energies are included in a RQRPA calculation, the substitution 

F(tt, t't', J n = + ) -> F(tt, t't', J n = 0+) Dm (45) 

is also pertinent. 

3 Results and Discussion 

The numerical calculations reported below are performed with two-body delta-force and 
the single particle energies used previously 0. We define the ratios 

v pp v pp 



pair ' pair ' 

Vs Vs 

between the T — 1, S — and T — (), S — 1 coupling constants in the pp channels and 
the pairing force constant, respectively. Their physical values are s — 1.0 and t = 1.5. An 
eleven dimensional model space, including all the single-particle orbitals from oscillator 
shells 3hu> and 4hu> plus O/2.9/2 and O/in/2 from the 5hu oscillator shell, is used. 
The results for two different RQRPA calculations will be discussed, namely: 

1. Approximation I (AI): the BCS equations are solved in the usual way with the 
constraint fl33|). 
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2. Approximation II (All): the pairing matrix elements are renormalized as indicated 
in ( fHf ) and the constraint (|54]) is used. 

As done in ref. [16], we analyzed the difference between the double iteration (DI) proce- 



dure, when the proton and neutron occupations are given by fllTl), and the simple iteration 
( SI) procedure, when the coefficients D pn on the right hand side of eq. (|TT| ) are put equal 
to unity. Besides, the influence of different number of multipolarities J n in the summa- 
tions involved in eqs. (|ll|) was studied by performing calculations either with only the 
J n = 1 + state (calculations denoted SI1 and DI1) or with the states J 71 " = l + ,2~,3 + 
and 4~ (calculations denoted SI4 and DI4)- The quantities that are of interest for the 
discussion are: the amount of ground state correlations, defined as 

C = Y,K + E^n, (46) 

p n 

the normalized GT strengths 

s _ = §_(r = i+)/(n - zy, s+ = s + (r = i + )/(n - z), (47) 

the normalized GT sum rule, 



S R = s_ - s+, (48) 

the lowest energy eigenvalue uj^ = i+ and the fifov matrix element given by ([21]) . 

As an example, we first discuss the calculations of these observables in 76 Ge nucleus. 
They are shown in Fig. 1 for different AFs (SI1, DI1, SI4 and DI4), as a function of 
the parameter t. The QRPA results, up to the collapse of this approximation, are also 
presented in the same figure. One can immediately see that: 

a) At variance with the QRPA, none of the RQRPA's collapses for physically mean- 
ingful values of t, and all four AFs yield quite similar results for u>jn = i+. 

b) The difference between the SI and DI procedures is always rather small. 

c) The GSC generated by the states with multipolarities J 71 " 7^ 1 + are not so significant 
as one would perhaps expect. 

d) The violation of the Ikeda sum rule is more pronounced when more GSC are taken 
into account. As the /3 + -strengths within the RQRPA and the QRPA are practically the 
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same, the deviation from the QRPA value S R = 1, mainly arises from the decrease of the 
/^"-strength. 

e) The different AFs yield very similar results for the (3f3 2v matrix element M.2v{J 1T = 
1 + ), and it is difficult to say which one is "better" and which one is "worse" . As a function 
of t, they all seem to behave in the same way as the corresponding QRPA matrix element. 
Still, the RQRPA results are qualitatively different in the sense that they cannot be fitted 
by the (1, 1)-Pade approximant 

M 2u ^M 2u (t = 0)^^, (49) 
i — t/ti 

as happens in the QRPA case ||. (t and ti denote, respectively, the zero and the pole 
of A4 2 u(J 7T = l + )0 This is a direct consequence of the absence of the collapse within the 
RQRPA for the physical values of t. 

All the above comments are also valid for the nuclei 82 Se, 100 Mo, 128 Te and 130 Te. 
From now on we discuss the differences between the AI and All results. In doing this 
we will recur only to the SI method. In Fig. 2 are compared the behaviors of the Ikeda 
reduced sum rule Sr for all five nuclei. | 

One immediately sees that in all the cases, except for 100 Mo, the violation of this sum 
rule is quite more pronounced in the AI than in the All. f\ Fig. 3 focuses the crossings 
of the calculated M. 2 v{J T: = 1 + ) matrix elements and the corresponding matrix elements 
deduced from the experimental half-lives (including the experimental errors) Q. Again, 
as we do not know which value of the parameter t should be used in each case, it is hard 
to say that the RQRPA is a better model than the QRPA. What is definitively clear is 
that the sensitivity of M.2v{J n = 1 + ) on this parameter is unavoidable and that it is not a 
consequence of the collapse of the QRPA. It has been pointed out more than once |19|, ^] 
that the Ai^^J 71 = 1 + ) amplitude goes to zero within the QRPA because of the partial 
restoration of the Wigner spin-isospin SU(4) symmetry. This is a result of two antagonistic 

2 It is worth noting that, when the SI procedure is used, Af n — J\f p and therefore S-(J n = + ) — 
S+( J 71 = + ) = (0|iV — Z|0). Thus from ( |32|) we see that the F sum rule is exactly conserved within All. 
Contrarily, in the AI it is broken by ^ 1% for s = 1 and by ^ 10% for s = 2. The F matrix elements 
M-2v{J 7! = + ) behave in the same way as the GT matrix elements and pass through zero at s = 1. 

3 From the theoretical point of view, the M.2v{J 7! = 1 + ) matrix element in 100 Mo is in same sense 
peculiar, because of the strong predominance of the [0^7/2 (n)0g9/ 2 (p); J" = 1 + ] configuration. 
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effects: the spin-orbit term in the mean- field potential that destroys the SU(4) symmetry 
(since it singles out one spin direction over the other) and the pn residual interaction that 
favors the LS coupling over the j — j coupling scheme [E3| . A physical criterion for fixing 
the triplet pp coupling strength t, based on the maximal restoration of these symmetries, 
has also been suggested (t sym ). Q In table 1 we compare the values of t sym , obtained using 
the recipe just mentioned, with those that are necessary to reproduce the experimental 
M-2v(J 7T = 1 + ) amplitudes (t exp ). The differences are of the order of 10% what is quite 
auspicious since the values of the nuclear parameter (both for schematic and realistic 
interactions) are not known with such a precision. We also remind that some additional 
degrees of freedom, not considered in the present calculations, such as the quadrupole 
and octupole charge- conserving vibrations can play an important role for all the nuclei 
discussed here. Besides, the contributions of odd-parity nuclear operators, arising from 
the p-wave Coulomb corrections to the electron wave functions and the recoil corrections 
to the nuclear current, are also significant for the /^/^-decays of 128 Te and 130 Te nuclei 
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In summary, when the GSC are taken into account the collapse of the QRPA does 



not develop in the physical region of the pp-strength parameter [16, O, 18]. Yet, the 



GSC only slightly mitigate the strong dependence of the j3p 2 u transition amplitude on 
this parameter, which is set on by the partial restoration of the SU(4) symmetry. The 
price that is paid to avoid the collapse within the QRPA is the non-conservation of the 
Ikeda sum rule within the RQRPA. This violation cannot be eluded and comes from the 
fact that the states generated by the action of the scattering part of the GT operator on 
the RQRPA ground state is not contained in the model space. 



4 By maximal restoration of the SU(4) symmetry we mean that for t = t sym the maximal concentration 
of the /3~-strength within the GT resonance takes place and the /3 + -strength is minimal. In no way it 
connotes that the J 71 " = 1 + states belong to a single SU(4) multiplet. 
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Figure Captions 



Figure 1: Behavior of several relevant quantities (defined in the text) within different RQRPA-yl/ 
approximations in the case the 76 Ge nucleus, as a function of t. The QRPA results, up to the collapse of 
this approximation, are also exhibited. The matrix elements M.iv are given in units of [MeF] -1 . 

Figure 2: Ikcda sum rule within the approximations 1 (upper panel) and II (lower panel), when the 
simple iteration method is used with only one intermediate state. 

Figure 3: Calculated matrix elements M.i v (in units of \MeV\~ 1 ) for several nuclei, as a function of 
t. Dashed and dotted curves correspond to the approximations I and II, respectively, and the solid lines 
are the QRPA results. The fill patterns indicate the experimental results. 



Table 1: Comparison between the values of the parameter t necessary to reproduce the 
experimental M.2u(J n = 1 + ) amplitudes (t exp ) and those that lead to maximal restoration 
of the SU(4) symmetry (t sym ). 
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